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ABSTRACT 



The probability of occurrence of extreme solar particle events (SPEs) with 
the fluence of (> 30 MeV) protons F30 > 10^° cm~^ is evaluated based on data of 
cosmogenic isotopes ^^C and ^^Be in terrestrial archives centennial-millennial time 
scales. Four potential candidates with F30 = (1 ^ 1.5) ■ 10-'^'^ cm~^ and no events 
with F30 > 2-10^" cm~^ are identified since 1400 AD in the annually resolved ^°Be 
data. A strong SPE related to the Carrington flare of 1859 AD is not supported 
by the data. For the last 11400 years, 19 SPE candidates with F30 = (1^3)- 10^° 
cm~^ are found and clearly no event with F30 > 5 • 10^° cm~^ (50-fold the SPE 
of 23- Feb- 1956) occurring. This values serve as an observational upper limit for 
the strength of SPE on the time scale of tens of millennia. Two events, ca. 780 
and 1460 AD, appear in different data series making them strong candidates to 
extreme SPEs. We built a distribution of the occurrence probability of extreme 
SPEs, providing a new strict observational constraint. Practical limits can be 
set as F30 ^ 1, 2-=-3, and 5 10^*^ cm~^ for the occurrence probability 10~^, 
10~^ and 10~^ year~^, respectively. Because of uncertainties, our results should 
be interpreted as a conservative upper limit of the SPE occurrence near Earth. 
The mean SEP flux is evaluated as fa 40 (cm^ sec)~^ in agreement with estimates 
from the lunar rocks. On average, extreme SPEs contribute about 10% to the 
total SEP fluence. 

Subject headings: Sun: particle emission — Sun: heliosphere — solar-terrestrial 
relations 
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Introduction 



Sporadic energy releases on the Sun can accelarate solar energetic particles (SEP) in 
the corona and interplanetary medium. Such phenomena often lead to solar particle events 



specifically Space Weather (iMewaldt 



(Mewaldt 


2006; 


Vainio et al. 


2009) 



20091 ) ■ A typical quantity for a 



SPE is the fluence of SEP with energy ab ove 30 MeV, Fsn- We have sufficient knowledge of 



20061 ) . with only several events with 



SPEs over the space era since mid-1950s (iSmart et al.l 
-^30 = (1 10) ■ 10^ cm~^ and hundreds of weaker SPEs observed. However, it is important 
to know, both for purely theoretical aspects of solar/stellar physics and for technical 
applications, the statistics of extreme SPEs with F-^q > 10^" cm~^. Such a study is possible 
only using indirect proxy data. So me estimates have been obtained from the measurements 



of cosmogenic isotopes lunar rocks (INishiizumi et al. 



20091 ) but it gives only the average fiux 



of SEP over very long scale without extracting individual SPEs. A list o f potential SPEs 



over the last 500 years was proposed base d on nitrate records in p olar ice (IMcCracken et al. 



2008 



2OI2I ). Cosmogenic nuclide ^^C 



2OOII ). but this result is heavily debated ( 1 Wolff et al. 
and ^'^Be data measured in terrestrial archives may provide informatio n on SPE in the past 



( ILingenfelter fc Hudson 



1980 



Usoskin et al. 



2006 



Webber et al. 



20071 ) but this possibility 



was not fully expl ored ear 



grossly uncertain ([Hudson 



ier. Accordingly, the probability of extreme SPEs remained 



2010h . 



Here we establish a solid observational constraint on the distribution of extreme SPEs 
using presently available datasets of cosmogenic isotopes (^"^C and ^"Be) measured in 
terrestrial archives with sufficient time resolution and quality, and modern models of their 
production in the atmosphere. We note that the presented resulted is based on terrestrial 
data and may not well represent the occurrence of solar events, whose geo-efficiency is also 
affected by the relative Sun-Earth attitude. 
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2. Data sets and the method 



2.1. Data 



The used data-sets are: 



a) IntCalOQ A^^C g lobal series: 11000 BC - 1900 AD, 5-yr time resolution 



( Reimer et al. 



20091). 



b) SB93 A^^C global annual series: 1511 - 1900 AD fIStuiver fc Braziunaslll993h . 



c) Dye3 ^°Be Greenland annual series: 1424-1985 AD flMcCracken et al.l 120041 ) 



d) NGRIP ^°Be Greenland annual series: 1389-1994 AD (IBerggren et alJl2009h . 



e) SP ^°Be South Pole Antarctic series: 850-1950 AD, quasi-decadal ( iRaisbeck et al 



1990 



Bard et al 



19971 ) 



f) DF Be Dome Fuji Antarctic series: 695-1880 AD, quasi-decadal (iHoriuchi et al 



20081 ). 



g) GRIP ^°Be Green land series: 7380 BC-1640 AD, quasi-decadal (lYiou et al.lll997 



Vonmoos et al. 



20061). 



h) M l 2 A^^C Japanese annual series: 750 - 820 AD, annual/biannual (IMiyake et al 
2012k 



2.2. Model computations 



In order to evaluate possible SPE signatures in the data, we used model computations 
of the isotope production by energetic particles, assuming instant injection of SEP into the 
atmosphere and calculating the expected isotope response in terrestrial archives. 



- 5 - 



As the referen ce event, we considered an extreme SPE of 23-Feb-1956 (SPE56 ) 



( iMeyer et al 



2009 



Usoskin et al. 



201l|). 



19561 ) with a very hard spectrum (iTylka &: Dietrich 
While SPE56 was the strongest observed ground- level enhancement (GLE), > 4000% in the 
count rate of the Leeds NM, it had relatively modest fiuence F^q = 10^ cm~^ (jShea fc Smart 



19901 ) ■ Yet sometimes SPEs with large fiuence but soft spectrum occur, e.g., a modest 
(only 10% at the polar Oulu NM) GLE event of 04-Aug-1972 (SPE72) with high fiuence 
F30 = 5 ■ 10^ cm~^. Since cosmogenic isotopes are produced by the most energetic part of 
the SEP energy spectra (> 1 GeV), we consider here the SPE56 scenario. 



Response of ^"Be to SPEs was ca 
using an updated ^°Be yield function (IKovaltsov fc Usoskin 



culated similar to lUsoskin et al.l (120061 ) but 



2OIOD, the corresponding 



geomagnetic model, and an intermediate at mospheric mixing 



Field et al. 



2006 



model (polar troposp heric 



Heikkila et al. 



200% . The 



and hemispherical stratospheric mixing - see 
calculated ^°Be production is 7.5 ■ 10^ atoms /cm^ for the SPE56 scenario. Because of 



the ^"Be stratospheric residence time (iHeikkila et al. 



20091 ). a possible SPE peak in ^^Be 
data can be 2-3-yr long. Since present models cannot convert the ^"Be production into 
conce ntrations in ice, they a re typically assumed to be directly proportional to each other 



le.g- 



McCracken et al. 



20041 ). For the normalization we use the reference period 1850-1900 



AD with moderate solar activity: th e galactic cosmic ray (GCR) modulatio n potential 



443 M V flAlanko-Huotari et al. 



20071 ) was close to the mean Holocene value (lUsoskin et al. 
20071 ). Accordingly, the reference period ^"Be production is (Q) = 3.45 ■ 10~^ at/cm^/sec 



or 1.1 ■ 10^ at/cm^/year. This value is scaled to the mean measured concentration for the 
reference period, in each ice core (1.025, 1.86, 3.47 and 2.76 10^ at/g for DyeS, NGRIP, SP 
and series, respectively). The GRIP series was normalized to the entire interval. 

The expected radiocarbon response of A^^C was calculated in two steps. First, the 



global production Q was computed using a new production model (IKovaltsov et al. 



2OI2I), 
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yielding 2.9 • 10 at/cm for SPE56. Next, the in stant atmospheric iniection of C (Q) 



2004h . and the 



was passed through a 5-box carbon cycle model (iDamon fc Peristykhl 
tropospheric A^^C was calculated. Fig. [T] shows the expected response in A^^C for SPE56, 
which is extended over decades and is characterized by a sharp increase and exponential 
decay, with the peak's FWHM 15-20 years. Accordingly, we look for a signature like that 
in the A^^C data. The SPE56 peak in the annual (5-year) A^^C data is 0.35 (0.19 ^ 0.27) 
permille, viz. below measurement errors of ~ 2 permille. Therefore, an SPE needs to be a 
factor XspE56 = 10 greater than SPE56 to be observable in the radiocarbon data. 

An SPE72-type soft event would require a 40-fold greater F^q, with respect to SPE56, to 
produce the same amount of cosmogenic isotopes. Accordingly, the estimates obtained with 
the reference SPE56 spectrum should be 40-fold enhanced to correspond to a soft-spectrum 
SPE72 scenario. 



3. Evaluating SEP fiuence in the past 

First we analyze data series with the annual resolution (b-d in Section [2?T!) . looking for 
potential signatures of SPEs. 

For the NGRIP series we looked for peaks with the magnitude > 1.3 ■ 10"^ at/g and 
duration < 3 years. Seven candidates were selected: 1436, 1460, 1650, 1719, 1810, 1816 and 
1965 AD. For the ^°Be Dye3 series five candidates (magnitude > 0.6 ■ 10'^ at/g, duration 
< 3 years) were selected: 1462, 1479, 1505, 1512 and 1603 AD. For the ^^C SB93 series 
no suitable candidates (peak with sharp rise and gradual decay) were found. Since the 
found candidates do not coincide in time, we performed a cross-check. For each candidate 
found in series A we calculated, using the model (Sect. 12.21) . the corresponding isotope 
production Q and the expected response in series B, and checked if it is consistent with the 
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data within ±2 years, and vise versa. For example, a 1436 AD peak in the NGRIP series 
must be accompanied by a strong peak in the DyeS series, which was not observed. Also, a 
NGRIP 1965 AD peak is rejected by direct cosmic ray observations. Using this cross-check 
seven candidates can be excluded: 1436, 1650 1816 and 1965 AD in NGRIP and 1479, 1512 
and 1603 AD in the DyeS series. The SB93 series is insensitive to check the ^"^Be-based 
candidates. Finally, five peaks pass through the check as potential SPE signatures (first 
block in Tabled]). Only one candidate (1460-1462 AD) is present in both series (discussed 
later) . 

Next we considered data with rougher time resolution. In each of the ^"Be data series 
(e-g in Section 12. ip with quasi-decadal resolution we searched for distinctive single peaks. 
Two candidates were identified in the GRIP series: ca. 5165 BC and 95 AD. One candidate 
was identified in the SP series ca. 1455 AD, probably related to the ca. 1460 event in 
NGRIP and DyeS annual series. Two candidates were identified in the DF series, ca. 780 
AD and 1805 AD. A search for signatures (sharp < 10-yr increase, FWHM 15-30 years, 
magnitude > 2 permille) in the IntCalOQ ^^C series yields many candidates (see an example 
in Fig. [3j\). We also made a cross-check of all the candidates between different series. For 
instance, the 1805 AD peak in the DF series is rejected basing on the annual NGRIP and 
DyeS data. About half of the candidates in the IntCal09 series after 7000 BC are rejected 
basing on the ^°Be data. Therefore, half of the identified candidates before 7000 BC can 
also be spurious. 

Finally, we list the candidates passing the cross-check, in Table [U One can see that 
a 10-fold SPE56 XspEse = 10 is the detection limit in cosmogenic isotope data. Only two 
candidates appear in more than one series - ca. 780 AD and 1460 AD, which are discussed 
below in detail. 
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Table 1: identified SPE candidates: approximate year, dataset, and tlie F30 fiuence [cm ^] 
evaluated for the SPE56 scenario. Scaling to SPE56 is given in parentheses. 



SPE year 


Series 


-^30 


(XsPEse) 


1460-1462 ADt 


NGRIP(1460) 


1.5 


• 10^° (15) 




Dye3 (1462) 


9.7 


■ 10^ (10) 


1505 AD 


Dye3 


1.3 


• 10^0 (13) 


1719 AD 


NGRIP 


1 • 10^° (10) 


1810 AD 


NGRIP 


1 • 10^° (10) 


8910 BC 


IntCal09 


2.0 


• 10^° (20) 


8155 BC 


IntCal09 


1.3 


• 10^0 (13) 


8085 BC 


IntCal09 


1.5 


• 10^0 (15) 


7930 BC 


IntCal09 


1.3 


■ lO^o (13) 


7570 BC 


IntCal09 


2.0 


• 10^0 (20) 


7455 BC 


IntCal09 


1.5 


• lO^o (15) 


6940 BC 


IntCal09 


1.1 


• 10^0 (11) 


6585 BC 


IntCal09 


1.7 


■ 10^0 (17) 


5835 BC 


IntCal09 


1.5 


• 10^0 (15) 


5165 BC 


GRIP 


2.4 


• lO^o (24) 


4680 BC 


IntCal09 


1.6 


■ W> (16) 


3260 BC 


IntCal09 


2.4 


• 10^0 (24) 


2615 BC 


IntCal09 


1.2 


• lO^o (12) 


2225 BC 


IntCal09 


1.2 


• 10^0 (12) 


1485 BC 


IntCal09 


2.0 


• 10^° (20) 


95 AD 


GRIP 


2.6 


• 10^0 (26) 


265 AD 


IntCal09 


2.0 


• lO^o (20) 


780 ADt 


IntCal09 


2.4 


■ lO^o (24) 




M12 


4 • 10^° (40)* 




DP 


4.5 


• 10^° (45)* 


1455 ADt 


SP 


7- 10^° (70)* 



Discussed separately. 
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50 

Years 



Fig. 1. — Response of the relative tropospheric concentration A^^C to SEP56 {t — shown 
by the arrow), with annual (grey curve) and 5-yr (black) time resolutions, depending on the 
5-yr bins boundaries. 



- 10 - 



3.1. Event of ca. 1460 AD 



Several analyzed series show a significant peak around 1460 AD (Tabl^I]). T 



his was 



noticed earlier and ascri 



2009 



McCracken et al 



3ed to a very strong SPE or t o a supernova explosion (e.g. 



I2004J : 



Delaygue fc Bard 



Berggren 



201ll ). We assume here that it is a SPE 



signature and evaluate its parameters. Annual NGRIP and DyeS series depict distinct 
peaks in 1460 and 1462 AD, respectively, which roughly agree with each other, and require 
-F30 ~ 10^° protons/cm^ for the SPE56 scenario. Rougher resolved SP series depicts a very 
strong peak ca. 1455 AD, requiring the fiuence 7 times greater than implied from the 
NGRIP/DyeS data. Neither of other series {GRIP, DF or IntGalOff) show peaks around that 
date allowing to evaluate the upper limit. Fig. shows a cross-check of this event vs. the 
DF series. While the event strength evaluated from NGRIP and DyeS series is consistent 
with no clear sig nal in the DF data, the huge fiuence implied from the SP data contradicts 



the DF data (cf. 



Delaygue fc Bard 



I2OIII ). A similar analysis of the GRIP data also suggests 
that the SP peak is too high. Fig. [2j3 shows the IntGalOO data along with signals expected 
from NGRIP/DyeS and SP data, respectively, for the ca. 1460 AD candidate. Because of 
the steep gradient in A^^C in 1430-1470, it is impossible to distinguish a small 2.4 permille 
peak implied by the DyeS/NGRIP data, but the high signal implied by the SP peak 
apparently contradicts the data. Thus, while all the analyzed datasets are consistent with a 
hypothesis of a strong SPE ca. 1460 AD, its fiuence is likely about 10^° protons/cm^ for the 



SPE56 scenario. The very high fiuence implied by the SP data contradicts other datasets. 



3.2. Event of ca. 780 AD 

Another interesting candidate for a strong SPE is ca. 780 AD as manifested through 
distinct peaks in DF, IntGalOO and M12 data series (Fig. [3]). According to these data-sets, 
the event was 25-50 times stronger than SPE56 with the fiuence F30 = (2^5)- 10^° 
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Fig. 2. — Test of a possible SPE ca. 1460 AD. A) DF ^°Be series, with points corresponding 
to the SEP-production estimated from NGRIP, DyeS and SP data series (see Legend) . B) 
IntCal09 A^^C series along with the expected signal estimated from DyeS/NGRIP and SP 
data series. 
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Fig. 3. — Test of a possible SPE ca. 780 AD. Open dots and black stars correspond to the 
modelled response to a SEP event equal to Xspe56 — 24 and XspEbs — 45, respectively. A) 
IntCalOg A^^C data series. B) M12 A^^C data series. C) DF i°Be data series. D) GRIP 
^°Be data series. 
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protons/cm^. The best fit for the IntCalOQ data (Fig. |3JA.) is obtained for a Xspe56 = 24 
event started ca. 780 AD. The M12 data (Fig. |3J3) is better fit with a Xspe56 = 40 event 
started in 775 AD, while the DF data (Fig. [SP) is consistent with Xspe56 = 45 event 
occurred between 780 and 790 AD. However, the GRIP series (Fig. [3p) depicts no peak 
at that period, which indicates that the fiuence greater than F^q = 3 ■ 10^° protons/cm^ is 
inconsistent at 0.03 significance leveL Applying the cross-check, we found that the fiuence 
implied by the IntCalOQ data {Xspe56 = 24 - see Tabled]) is consistent with the GRIP 
data, while both Xspe56 = 40 and 45 based on the DF and M12 data, respectively, yield a 
too high peak in the GRIP series. Accordingly, we conclude that the event of 780 AD is a 
possible candidate for a strong SPE with a consensus value F^o ~ 3 ■ 10^° protons/cm^. 



3.3. Carrington event of 1859 AD 



The Carrington event of 1859 AD is often considered as an extreme SPE, with 
■ F-^n = 1.8 • 10^° cm ~^ as estimated from the nitrate record in the Greenland Summit core 



(IDreschhoff Sz Zeller 



1998 



McCracken et al. 



2001 



Shea et al. 



20061 ). Although none of the 



series analyzed here depicts a peak around 1859, we check if this prop osed SPE is consistent 



with the cosmogenic data. Using the G CR modulation reconstructio n (lAlanko-Huotari et al. 



20071 ) and the ^"Be production model (iKovaltsov fc Usoskin 

expected ^"Be concentration along with the additional production for 1859 AD, app 



20101 ) ■ we calculated the 



McCracken et al. 



(1200lh . 



SEP parameters for the Carrington SPE according to Section 8 of 
The model computation is then confronted with the annual data of two Greenland series, 
NGRIP and DyeS, in Fig. HI The ^''Be peak, expected from the nitrate data, is too strong, 
contradicting to the observed data in two Greenland sites. Radiocarbon data and decadal 
^"Be series cannot resolve the Carrington peak. Therefore, we conclude that the cosmogenic 



data do not support the hypothesis of a very strong SPE related to the Carrington fiare. 
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Fig. 4. — Testing a possible SEP event ca. 1859 AD. The gray line depicts a model signal 
by GCR plu s an addition of the ann ual production due to the Carrignton 1859 AD event 

(1200 li l. A) Time series of the annual NGRIP ^°Be data, with 



according to 



McCracken et al. 



15% error bars. B) Time series of the annual DyeS ^^Be data, with 15% error bars. 
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4. Discussion 

We propose, by analyzing data from different cosmogenic isotope records and 
performing tlieir cross-check, a list of candidates for strong SPEs in the past (Tabled]). The 
list is somewhat uncertain, as some events may be missing or some spurious peaks falsely 
identified. However, the result is quite consistent in a statistical sense. 

Four candidates with the Xspe56 factor 10-^15 are identified in the annual ^°Be series 
NGRIP and DyeS for the last 600 years (first block of Table [T]). Events with XspEse < 10 
cannot be reliably identified. We can also securely say that events with XspEse > 20 are 
not observed during that period (cf. Fig. H]). In particular, a strong SPE related to the 
Carrington flare in 1859 contradicts these data. 

We evaluate the probability p of occurrence of such events, applying the Poisson 
distribution (assuming that SPEs are mutually independent). For example, for four events 
observed during 600 years, the probability is p = 0.00771o;oqJ5 year~^ (90% confidence 
interval), which is higher than the naively taken 4/600 = 0.0067 year"^. If no event with 
-F30 > 2 ■ 10^*^ cm~^ is observed over 600 years, the corresponding one-sided 90% confidence 
interval is p = ^ 0.0027 year~^, with the median probability 0.0012 year~^ (once per 850 
years) . 

The last 4 -6 centuries co ver the full range of solar activity, from Grand minima to 



the maximum ( lUsoskin 



20081 ) serving as an archetype of the solar variability. Therefore 
we expect that this result is consistent with longer scales. On the multi-millennial scale 
(the Holocene - 11 millennia), 19 candidates are identified (second block in Table [1]) with 
-^SPE56 = 10 ^ 30. A few candidates with XgpEse > 50 were rejected by the cross-check. 
This gives the average occurrence rate of such SPEs roughly two per millennium. Moreover, 
we can securely say that during the Holocene there was no events with XgpEse > 50 100, 
placing a strong upper limit on the SPE strength. We found no apparent relation of the 
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SPE candidates occurrence and the solar activity level (lUsoskin et al.l 120071 ) 



We summarize our findings in a plot (Fig. E]) of the integral probability of strong 
SPE occurrence. X- and Y-axis give the F^o fluence and the probability of occurrence 
of SPE with the > 30 MeV proton annual fluence exceeding F^q, for different datasets. 
he measured annua l fiuences for the space era (1956-2008) are shown by open triangles 



(IShea fc Smart 



1990l . Shea, pri vate communication 2012). Since the space era coincides 



with the unusually active Sun f ISolanki et al. 



20041 ). this probabilities may be higher than 



for the typical medium activity Sun. The black triangle reflects the fact that no SPE with 
F30 > 10^" cm~^ was observed during 53 years. The open star corresponds to the four 
SPE candidates from the annual ^"Be data (first block in Table [1]), while the filled star 
corresponds to no SPE with F^q > 2 • 10^° cm~^ found during 600 years. Open circles 
represent the candidates found in the rougher time series (second block in Table [I]), grouped 
into four points F^q > 1, 1.5, 2 and 2.5 10^*^ cm~^, respectively. Since the lower fluence 
value is probably underestimated because of the detection threshold, we ignore the first 
open circle, using instead the star-symbolled point from the annual ^°Be datasets. The 
filled circle corresponds to no event with F^q > 5 • 10^° cm~^ found over 11 millennia. One 
can see that there is a break in the di stribution at about -F30 = 5 ■ 10^ cm~^. This break (c f. 



McCracken et al 



2001 



Hudson 



2OIOI ) may be related to a streaming linait ( 



and/or limitations in the energetics of solar flares/CME (IFletcher et al 



R 



eames 



19991) 



20111 ). Our new 



results put a stronger constraint on the distribution. For example, a SPE XspEse > 20 is 
expected only once per millennium, and no XspEse > 50 event is expected over tens of 
millennia. We try to fit the tail of the distribution (F30 > 7 • 10^ cm~^) by two simple 
models: exponential (p oc exp (—0.33 ■ 10~^ ■ -F30)) and power law (p oc -^30^) - see grey lines 
in the Figure. The data does not allow us to distinguish between the two shapes. 



We also compare our results with the average SEP flux over time scales of thousands- 
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Annual fluence F [10^ cm"^] 



Fig. 5. — Probability of occurrence of the annual fluence (> 30 MeV) exceeding the given 
value F30, as evaluated from the data for the space era 1956-2008 (triangles), cosmogenic 
isotope annual data (stars) and cosmogenic isotope decadal data (circles). Open symbols 
correspond to the statistics of "observed" data, while filled symbols denote the conservative 
upper limit (no events greater than X observed during Y years) . Error bars correspond to the 
90% confidence level evaluated assuming the Poisson statistics. Grey dotted and solid lines 
depict the best- fit exponent and power law approximations for F30 > 7 • 10^ protons/cm^. 
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millions years, estimated fro m lunar rocks, which ranges 21^56 (cm^ sec) ^ (e.g. 



1999 



Nishiizumi et al. 



Reedy 



20091 ). Time averaging of the results from Fig. [S] using the above fits 
for the tail yields the average flux ^ 38 (cm^ sec)~^, whi c h is compose d of about 35 (cm^ 



sec) ^ for the space era record (e.g., 



Shea fc Smart 



2002 



Reedy 



2OI2I ) with an addition 



of 3.2 (cm^ sec)~^ due to the distribution tail based on cosmogenic isotope data (both 
exponent and power law give similar results), suggesting that extreme SPEs contribute 
only about 10% to the total SEP fiuence. This is totally consistent with the assessments 
based on lunar rocks, giving an independent support to the present results. For a SPE 
with a softer energy spectrum, like SPE72, the tail addition would be 40-times greater 
(Section 12.21) leading to the average SEP flux ^^^150 (cm^ sec)~^, in contradiction with the 
lunar rock data. Therefore, extreme SPEs found in the cosmogenic isotope records must 
have (on average) hard energy spectra. 



5. Conclusions 

We have evaluated the probability of occurrence of extreme SPEs based on data of 
cosmogenic isotopes ^^C and ^"^Be in terrestrial archives, spanning over the time scale from 
centuries to 11 millennia (FigE]). We identified four potential candidates for SPEs with 
F30 = (1 ^ 1.5) ■ 10^° cm~^, and show that no event with F30 > 2 ■ 10^" cm~^ existed, over 
the last 600 years using annually resolved ^"Be data. In particular, the extreme Carrington 
SPE of 1859 AD contradicts these data. From rougher resolved data we identified 19 SPE 
candidates (Tabled]) with F30 = (1^3)- 10^*^ cm~^, and clearly no event with F30 > 5 ■ 10^" 
cm~^, over the last 11400 years. Two events appear in different series, ca. 780 AD and 1460 
AD making them strong candidates to extreme SPEs. This gives a new strict observational 
constraint on the occurrence probability of extreme SPEs. Practical limits can be set as 
F30 ^ 1, 2 ^ 3 and 5 10^° cm'^ (10-, 20^30- and 50-times greater than SPE56), for the 
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occurrence probability of 10~^, 10~^ and 10~^ year~^, respectively. The mean SEP flux is 
found as ^ 40 (cm^ sec)~^ in agreement with estimates from the lunar rocks. On average, 
extreme SPEs contribute about 10% to the total SEP fluence. 

We note that the present result tends to represent an upper hmit for the SPE 
occurrence, since we explicitly assume that every peak in one data series, consistent 
with other series, is a signature of SPE. This may be not exactly correct as some of the 
peaks may be spurious. Accordingly, our results should be interpreted as a conservative 
upper limit of the SPE occurrence near Earth. Application of the results to the Sun is 
not straightforward, especially for most energetic events with low statistics, because the 
propagation of SEP in the interplanetary space may greatly affect geo-efficiency of SPE, 
and accordingly, their ability to become detectable in the cosmogenic isotope data series. 
Given these uncertainties, the present results should be considered with a precision of up to 
an order of magnitude. 
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supported by the Program No. 22 of the Presidium RAS, and acknowledges also support 
from the Academy of Finland and Suomalainen Tiedeakatemia ( Vaisala foundation) . 
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